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4 Abstract. 

5 The magnetopause boundary implies two kinds of variations: a density/ 

e temperature gradient and a magnetic field rotation. These two kinds are al- 

7 ways observed in a close vicinity of each other, if not inseparably mixed. We 

s present a case study from the Cluster data where the two are clearly sep- 

9 arated and investigate the natures of both layers. We evidence that the first 

io one is a slow shock while the second is a rotational discontinuity. The inter- 

ii action between these two kinds of discontinuities is then studied with the help 

12 of 1,5-D magnetohydrodynamics simulations. The comparison with the data 

13 is quite positive and leads to think that most of the generic properties of the 
« magnetopause may be interpreted in this sense. 
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1. Introduction 

is The Earth magnetopause is the outer boundary of the terrestrial magnetosphere. Out- 

i6 side of this boundary, the magnetosheath plasma is the shocked solar wind plasma, i.e. 

i? cold and dense, with a magnetic field direction essentially determined by the solar wind 

is one. Inside of it, the magnetospheric plasma is hot and tenuous, with a magnetic field 

lo direction essentially determined by the planetary one. For this reason, two kinds of strong 

20 gradients are observed in the same close vicinity: a compressional one (density, temper- 

2i ature and possibly magnetic and kinetic pressures) and a rotational one (magnetic field 

22 direction). Understanding how these two kinds of variations can co-exist is a pivotal issue 

23 for understanding the nature of the boundary. Whenever this boundary can be locally con- 

24 sidered as stationary and one-dimensional, the Rankine-Hugoniot jump equations allow 

25 only one type of discontinuity able to ensure simultaneously the two kinds of variations: 

26 the tangential discontinuity. It implies that the normal magnetic field and flow are equal 

27 to zero in the frame of the discontinuity. 

28 There is now observational evidence (Chou et al, 2012) that the magnetopause is not al- 

29 ways a tangential discontinuity, i.e. that the normal flow and magnetic field are not always 

30 null. There are in particular numerous observations where the magnetic field rotation cor- 

31 responds to a rotational discontinuity (Sonnerup et al, 1974). In this paper, we use the 

32 word "discontinuity" in its usual sense in this context: a 1-D stationary layer, whatever its 

33 thickness. It seems that, at least when dynamical processes takes place (due for instance 

34 to reconnection near the studied site), the compressional boundary and the rotational one 

35 can be distinguished, both propagating with different velocities with respect to the flow. 
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36 Our goal is to present evidence of such a case with a detailed experimental analysis and 

37 suggest an interpretation of the observations. In sections 2 to 4, the experimental Clus- 

38 ter case is presented and the two crossings are analyzed with the help of a new method 

39 presented in [Dorville et al, submitted to JGR 2013]. In section 5, the interaction of a 

40 compressional jump and a rotational one is investigated by solving numerically the MHD 
4i equations, in order to understand to what extend the properties of both are modified by 

42 the interaction, in particular their propagation speeds. The comparison with the data 

43 and the conclusions are drawn in section 6. 

2. Study of a magnetopause crossing by Cluster 

44 Using the single-spacecraft method of analysis described in (Dorville et al, submitted 

45 to JGR, 2013) we are able to build a one-dimensional spatial transition parameter at the 

46 magnetopause and to analyze the small scales inside the boundary during a magnetopause 

47 crossing. This new transition parameter is expectingly close to a normal linear coordi- 

48 nate. The assumptions of the method is that the boundary is one-dimensional, sufficiently 

49 stationary during the crossing, and that the speed of the flow in its frame is negligible 
so with respect to the speed of the boundary. A fit of the two temporal tangential field 
si components with a spatial elliptic form is done, in the frame where the normal magnetic 

52 field is constant and where the normal speed of the flow is consistent with an angle on 

53 the ellipse linear with spatial position. 

54 

55 In this part we focus on a magnetopause crossing by Cluster CI on April 15th 2008, 

56 between 15:19:05 and 15:25:05. Fig. 1 shows the variation of ion density obtained by the 
s? Hot Ion Analyzer (HIA), the energy spectrum, and the variation of the Fluxgate Magne- 
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58 tometer (FGM) magnetic field (Balogh et al, 1997) associated with this crossing. 

59 For this magnetopause crossing, we are able to fit the magnetic field and find a transition 

60 parameter consistent with CIS (Reme et al, 1997) ion velocity. The observed angle on the 
6i ellipse as a function of time (in black) and the angle derived by the fit with CIS velocity 

62 (in red) are plotted on Fig. 2. The fit is good on this case. The method also provides 

63 an accurate determination of the normal direction, which is here slightly different from 

64 the MVABC [Minimum and Maximum Variance Analysis] one: ([0.98,-0.18,-0.10] against 

65 [0.96, -0.10, 0.25] in GSE coordinates). The result of the fit on magnetic field is presented 

66 on Fig. 3, where time zero is the center of the large data interval and the fit is done on the 
e? time period of the magnetic field rotation. It is worth noticing that the variance on B m 

68 is not much larger than the variance on B n , which makes the MVAB method inefficient. 

69 On the other hand, MVABC assumes that B n = 0, which does not allow determining it 

70 either. Our method provides a value B n = 14nT, corresponding to a propagation angle 

71 @Bn = 75 degrees . 

72 

73 Looking inside the boundary, Fig. 4 presents the time evolution of density, total magnetic 

74 field and Bz. It seems clear that there are two phases in the evolution of these quantities: 

75 first, the density is rapidly increasing and the magnetic field amplitude B is decreasing. 

76 Then, these quantities remains constant and the magnetic field slowly rotates on an el- 

77 lipse. We will now study these two phases separately, assuming they are oriented along 

78 the same normal direction. 
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3. Compressional variations 

79 Let us first focus on the first, compressional, phase, between t=-9 and t=2, when the 

so density is growing and the total field decreasing. The first step to study this discontinu- 

8i ity is to find its normal velocity, in order to work in its proper frame. We cannot here 

82 make the simple assumption that the observed velocity is the boundary velocity because 

83 this layer is compressional and there is certainly, according with the conservation of the 

84 flux, a consistent velocity gradient inside the boundary. To determine the velocity of the 

85 boundary using the mass conservation, we plot the measured ion flux pu n relatively to 

86 the measured density p Fig. 5, for the three CIS measurement points available inside the 
a? shock, altogether with the nearest points around the shock. This flux has to be constant 
88 in the frame on the boundary, and its variation must be proportional to the velocity of 
as the boundary in another frame. We find that the points are aligned with a constant 
90 flux of 1.07 ± 0.11 10 11 (SI) and a velocity of the boundary of 55.1 ± 1 Km/s, the given 
9i uncertainties being the 1-sigma uncertainties estimated from the fits. If we restrict to the 

92 only three points in the gradient, we find that the points are perfectly aligned with a flux 

93 of 1.09 ± 0.01 10 11 (SI) and a velocity of the boundary of 54.7 ± 0.1 km/s. 

94 

95 Here we have to notice that if we look at the same points on MVABC frame or MVAB 

96 frame, the flux is not proportional to the density. So the precision of the frame found with 

97 our method seems to be better, that is very important for such a precise treatment. Even 

98 with the small numbers of CIS points, we can be confident on this value and work on the 

99 discontinuity frame. Fig. 6 represents the variation of density, pressure, field modulus 
wo and normal velocity. We see that the last two quantities decrease as the others increase, 
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that is possible only for a slow shock with respect to the Rankine-Hugoniot equations. 

Indeed, we find that the flow is super-slow on magnetospheric side and sub-slow on the 

other side, as the Cs on the curve refers to the slow mode velocity. So this discontinuity 

can be identified with confidence as a slow shock. 

4. Rotational variations 

Fig. 7 shows the second part of this magnetopause crossing. The plotted quantities are 
the density, temperature, pressure and magnetic field. All seem to be constant in this 
part of the boundary, with only a slow rotation of the direction of magnetic field. It is 
also the case for the velocity modulus. These evolutions are characteristic of a rotational 
discontinuity. 
To test this discontinuity, we used a variant of the Walen test (Paschmann et al, 2008) 



in on Fig. 8. We plot the variations 5B t as a function of J(noP) $ u t for the two tangential 

U2 components. The result appears consistent with a rotational discontinuity, but with a 

us little difference: the slope is not equal to one as it should, but to about 0.72, 0.74 with 

U4 the correction coming from the temperature anisotropy at this place. To explain it from 

us the composition of the plasma, we would need a mean mass of 1.8 protons. Unfortunately 

U6 there is no CODIF data for CI on April 15th, and no good resolution electron density 

a? measurement. 
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5. Direct numerical simulations of a slow-shock/rotational discontinuity 

interaction 

120 Let us now consider the hypothesis that this magnetopause crossing is the observation of 

i2i a non stationary interaction between an Alfven wave and a slow shock. This model should 

122 enable to describe many observations of the magnetopause, since both a compression and 

123 a rotation of B have to be present. In order to investigate this possibility, we use a 

124 one-dimensional MHD simulation code. 

5.1. Initialisation of the shock 

125 We use a 1.5 D compressible MHD simulation code with periodic boundary conditions. 

126 The code uses a classical Fourier pseudo-spectral method to compute spatial derivatives, 

127 with an adaptative second order Adams-Bashforth time stepping for time integration. 

128 We integrate the equations for the density p, the gas pressure P, the velocity u, and the 

129 magnetic field B. The temperature is defined by an ideal gas law. These equations read, 

130 appropriately normalized: 

d t p + V-(pu) = (1) 

d t P + (u ■ V)P + 7?V-m = Di (2) 

d t u + (u ■ V)u + V(P + B 2 /2)/p - (B ■ V)B/p = D 2 (3) 

d t B + (u ■ V)B -(B-V)u + BV-u = D 3 (4) 

P = pT (5) 

i3i D 2 and D 3 denote respectively the viscous and resistive dissipative terms, and D\ the 

132 corresponding heating terms, together with conduction. Only plane waves are considered, 

133 i.e., V = (d x ,0,0). The first step for this study is to obtain a sufficiently stable slow 

DRAFT April 5, 2013, 1:34am DRAFT 



DORVILLE ET AL.: MAGNETOPAUSE: ROTATIONAL/COMPRESSIONAL X - 9 

134 shock, in order to be able to see the interaction with an Alfven wave. 

135 To determine the needed jumps between the two sides of the shock for all parameters, we 
we use the Rankine Hugoniot equations. We computed the different possible couples u n i,u n2 
137 for the upstream and downstream normal velocities (see Fig. 9), for any values of the in- 
ns cident magnetic field and pressure. The curve presented on Fig. 9 is for example obtained 
i3, for Pi = 0.1, B n = 0.1, pi = 1, and B tl = 1. 

«o We inject the chosen values u n \ = 0.065, u n 2 = 0.035 in the simulation code, as the 

i4i extreme values of a ramp of u n with a tanh kx profile. The width of the shock is of a 

«2 few mesh points. To respect the periodic boundary conditions, we add a linear variation 

143 leading the u n function to the same value on the two boundaries. To initialize all the other 

144 quantities at all points, we directly use the conservation laws of the Rankine-Hugoniot 

145 system. 

«6 Using this initialization we obtain something which is not perfectly stationary (our equa- 

147 tions do not take account of the code viscosity and the profiles are thus not perfectly 

us known), but enough to study the interaction with a rotational discontinuity. Fig. 9 shows 

us the evolution of density and fields for small times. A small perturbation of density seems 

150 to propagate alone on the right on the box (we will see that it does not change the results), 

i5i and we observe a small smoothing of the gradients, but the shock is not globally moving. 

5.2. Interaction with an Alfven wave 

152 We then add an Alfven wave in the right side of the simulation box, consisting in a 2ir 

153 rotation of the magnetic field , with a Walen consistent rotation of velocity. The subse- 

154 quent evolution after this initialization is shown on Fig. 11 for small times (until t=15). 

155 We see that the Alfven wave is normally propagating on the x direction from right to left 
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and encounters the density gradient around t=9. The shock then begins to move with the 
Alfven wave instead of remaining steady. The imperfection on the initialization does not 
play any significant role in this phenomenon. 

This effect becomes spectacular if we look at the simulation for a longer time. Fig. 12 
and Fig. 13 show the evolution every 20 seconds during 60 seconds, respectively with the 
Alfven wave, and without the Alfven wave. There is no doubt that we see the formation of 
a slow shock/rotational discontinuity compound, the slow shock moving with the Alfven 
wave. Such a structure is quite reminiscent of the so-called "double discontinuities" de- 
scribed in (Whang et al, 1998). Nevertheless, these structures are usually supposed to be 
related to an anisotropic (Lee et al, 2000) plasma, which is not the case in this MHD code. 
However, the anisotropy is said necessary to explain the stationarity of this structure. As 
the observed layer is not strictly stationary, even if it evolves very slowly when the two 
layers merge, there is no necessary contradiction. 

This kind of structure is forbidden by MHD equations if it has to be stationary, but it 
is not, at least at times we can reach with the simulation, as shown on Fig. 14 that shows 
the respective positions of the Bz and density gradients at different times. The shock 
affects the shape of the rotation of magnetic field (even if a Walen test can show it is still 
Alfvenic), and the compressional part position seems to be oscillating. 



6. Conclusion 

We have presented a case study of a magnetopause crossing by Cluster spacecraft. 
This example is atypical in the sense that the density/ temperature gradient is clearly 
separated from the magnetic field rotation layer. This allows to investigate separately 

DRAFT April 5, 2013, 1:34am DRAFT 



DORVILLE ET AL.: MAGNETOPAUSE: ROTATIONAL/COMPRESSIONAL X - 11 

its the two layers and determine their nature. Using a new method, the normal direction is 

179 determined with a good accuracy and the normal components B n of the magnetic field 

wo and u n of the velocity can be measured in a reliable way: the two layers are not tangential. 

isi In addition, the method allows plotting the different profiles in function of a transition 

182 parameter which is a trustworthy proxy of a normal coordinate. The density /temperature 

183 gradient layer could so be identified unambiguously as a slow shock, while the rotation 

184 layer has all the properties of a rotational discontinuity, except that it may propagate at 
las 0.7 Va instead of Va- A 1,5 D MHD simulation has been performed to investigate the 
we interaction between a slow shock and a rotational discontinuity when they meet. The 
is? results are quite comparable with the data : The shock stays attached to the rotational 

188 discontinuity, at least before viscosity kills all the gradients in the simulation, and keeps 

189 shock properties as the second part of the discontinuity is still rotational. The structure 

190 is not stationary and so stays compatible with MHD equations. We can conjecture that 
mi this kind of interaction between compressional and rotational features is a generic feature 

192 of the magnetopause. 

193 Acknowledgments. The authors would like to thank the CAA and all Cluster instru- 
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Figure 1. Density, energy spectrogram and magnetic field observed by Cluster CI around 
15:21 on 04/15/2008. The jump of density, change in plasma energy composition and rotation 
of magnetic field show that the satellite is crossing the magnetopause. 




time(s) 
Figure 2. Observed angle on the elliptic hodogram and fit (dashed line) derived with CIS 

velocity by the BV method as a function of time for Cluster CI crossing of 04/15/2008. 
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Figure 3. Three components of the magnetic field data measured by FGM and fit (dashed line) 
by the BV method in the magnetic field rotation region for Cluster CI crossing of 04/15/2008. 
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Figure 4. Time evolution of density, total magnetic field and B^ component of the magnetic 
field during the for Cluster CI magnetopause crossing of 04/15/2008. 
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Figure 5. The normal mass flux in the shock (stars) and around it (crosses) for the BV method 

improved frame and in the MVABC frame. The proportionality in the BV frame shows that it 

is considerably better and permits to deduce the shock velocity with respect to the spacecraft. 
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Figure 6. Normal velocity, density, pression, magnetic field modulus and slow mode velocity 
as a function of time for the shock. 
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Figure 7. Bl component of magnetic field, density, temperature, and magnetic field modulus 
as a function of time for the rotational discontinuity. 
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Figure 8. Variant of the Walen test for the rotational discontinuity, y- * dVt is represented 
over dB t . 
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Figure 9. Possible couples of normal velocity boundary conditions allowed by Rankine- 
Hugoniot equations (Vn2 as a function of Vnl). 
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Figure 10. Early evolution of plasma parameters from the initialization (blue curve) to the 
last time step (red curve) for the shock initialized alone in the box. 
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Figure 11. Early evolution of plasma parameters from the initialization (blue curve) to the 
last time step (red curve) when a rotational discontinuity is added in the initial condition next 
to the shock. 
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Figure 12. Long time evolution of plasma parameters from the initialization (blue curve) to 
the last time step (red curve) when a rotational discontinuity is added in the initial condition 
next to the shock. 
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Figure 13. Long time evolution of plasma parameters from the initialization (blue curve) to 
the last time step (red curve) for the shock initialized alone in the box. 





Figure 14. Shock and rotation discontinuity positions at different times. The density is used 
as a proxy for the shock and the Bz component as a proxy for the rotational discontinuity. The 
structure is not stationary. 
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